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Abstract 

The Feynman-Schwinger representation is used to study the behavior of so- 
lutions of scalar QED in (2+1) dimensions. The limit of zero photon mass is 
seen to be smooth. The Bethe-Salpeter equation in the ladder approximation 
also exhibits this property. They clearly deviate from the behavior in the non- 
relativistic limit. In a variational analysis we show that this difference can be 
attributed to retardation effects of relativistic origin. 
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An important issue in the description of relativistic composite systems is to 
find a reliable and practical formalism that is consistent with known limiting 
cases and the underlying symmetries of the system at hand. The Feynman- 
Schwinger representation (FSR) [1-7] offers a nonperturbative description 
which is consistent with both gauge and Lorentz invariance. Furthermore, the 
formalism was shown to satisfy the correct static and nonrelativistic limits. In 
this formulation all ladder and crossed graphs are summed up, while the in- 
clusion of valence particle loops, self energies and vertex corrections is feasible 
as well. Moreover, the FSR is well suited for the study of nonabelian confining 
theories such as QCD, since vacuum condensates can easily be accounted for 
via the cumulant expansion [3-5]. These considerations indicate that the FSR 
is an appealing alternative to the celebrated Bethe-Salpeter equation (BSE). 

In this letter we apply the FSR to (2+1) dimensional scalar QED (SQED2+1) 
in the two-body sector. Our interest in this theory is twofold. First we consider 
it to be an excellent playground to study the FSR in the equal mass case, since 
the theory is UV finite in the two-body sector in contrast with the situation 
in (3+1) dimensions. Secondly, it has been proven that compactified quenched 
spinor QED 2 +i is linear confining for all values of the charge ||[| . Since the 
spin of the valence particles is usually not considered to be essential for the 
mechanism behind confinement, one may hope to learn something about the 
confinement mechanism in this abelian case. In this work however, we will 
focus on the restoration of gauge invariance. 

In order to avoid IR problems we introduce a photon mass \x in the theory 
and study the limit — > 0. In this limit we may compare the FSR-results with 
those of the BSE in the ladder approximation and the Schrodinger equation. 
The \i dependence can readily be found in the latter case. Taking the one 
meson exchange contribution as driving force between 2 charges with mass m, 
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we have in 2 spatial dimensions: 



1 e 2 \ 
—A - —K (fj,r) ^(0 = -E B ^(r) 
m 2tc I 



As \i — » the modified Bessel function Ka(/j,r) behaves as — log(/xr), which 
leads to the well known result of a logarithmic confining Coulomb potential 
in 2 spatial dimensions after an (infinite) energy shift. This interaction causes 
E B to diverge as — |^log(/i/ey / m), at least for the low lying states. 

Let us now turn to the sQED 2 +i case. We consider 2 scalar particles with mass 
m, minimally coupled to the massive photon field A^. The Euclidean action 
for this theory is: 



S 



J d 3 x [|(<9, - ie^)0| 2 + m 2 |0| 2 + \Fl + \^A\ + \C\d,A,f] (2) 



The parameter £ takes care of the gauge fixing when we restore gauge invari- 
ance by taking the limit fi — > 0. 

The object under study is the gauge invariant 4-point function of the theory, 
defined as the transition matrix element between the initial state ^f?i(x, x) = 
(fy(x)P(x, x)4>(x) and the final state ^f(y,y) = (f)(y)P(y,y)(p i (y): 

G(x, x, y,y)= J V<P VA^ * f (y, y) %(x, x) e~ s (3) 



The wavefunctions and \l/f are defined in a gauge invariant fashion by means 
of the parallel transporter P(x,y) = exp[— ie f% dz^A^z)). Next the valence 
field (f) is integrated out and the resulting determinant is set to unity. The 
latter amounts to neglecting all 0-loops. It was shown in |I| that the resulting 
'quenched' Greens' function can be written in the following form: 

oo oo 

G(x, x, y,y) = JdsJdsJ {Vz) xy {Vz) m e -*tM,*] (w {Ktg} ) (4) 





3 



where the path integral / (T>z) xy is a quantummechanical one, subject to the 
condition of fixed endpoints: z(0) = x and z(l) = y. The functional K is the 



action of a free relativistic point particle with fixed eigentime s |TT 

l r 

K[{z},s]=m 2 s + -j drij(r) 



(5) 



The Wilson-loop (W{z,s}/ can be computed exactly in the [/(l)-case consid- 
ered here: 



W{*,*}> = ( exp 



exp 



-ie y dw M A M (w) 
c 



-~e j> dv^ j dwuA^v - w) 

c c 



(6) 
(7) 



W{ z ^z} ) is known to be an order parameter of the deconfinement transition in 
pure gauge theories jlO|JTTi . The variables v and w in @ and (0) both run along 
the whole closed contour C formed by the lines xx and yy and the paths {z} 
and {z}. Any longitudinal term vanishes identically in (0)0- The nonvanishing 
transverse component of A^ u is given by: A AtI/ (r) = 5 /Jiy e~' ir '/47rr = 5 MI/ A(r). It 
can be shown that since there is no ordering of t> and w present in (|7|), it sums 
up all ladder, crossed ladder, self energy and vertex correction graphs with 
the appropriate weights. Comparing the FSR to quenched lattice gauge cal- 
culations, the calculation of (F§) can obviously be carried out in a numerically 



3 In order to see this we Fourier transform to momentum space and concentrate on 
the longitudinal contributions: 



dr / dcr w^(t) v v (a) A^ iu (w - v) 
d d q 



dr / da d T (q ■ w(r)) d a (q ■ v(a)) f(q) 



(2tt) c 

which vanishes identically for every closed contour. Notice that this is the case for 
finite jJL as well. 
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more accurate way since the occurring path integrals are quantummechanical 
ones. 



Here we wish to concentrate on the ladder and crossed graphs only and we 
therefore restrict w and v to opposite sides of the contour. The bound state 
spectrum of the thus obtained expression is studied by considering the asymp- 
totic spectral decomposition of G: 

G(T) = lim £ c n exp(-£ n T) (T = ±(y 3 + y 3 - x 3 - x 3 )) (8) 

1 — »oo 

n 

and its logarithmic derivative L(T) = - [^G(T)1 /G(T). According to (§), 
L(T) is expected to approach the mass of the ground state Eq as T —>■ oo. In 
practice we approximated the path integral (|J) by a finite product of iV — 1 
integrals: 

r / N \ 3N/2 r 

with the constraints z = x and Zfj = y. The normalization factor in (^) is not 
irrelevant since it contains s which is integrated over. The functionals K and 
V = log (W) are discretized as follows: 



N N 

K[{z}, S ]^m 2 s + -Y, & ~ z t -if (10) 
4s i=i 



V[{z,z}} 



N 

e 2 



( z i - Zi-i) ■ ( z j - Zj-i) A (| (zj + zt-i - Zj - Zj-x)^ (11) 



;j=i 



Here we adopted the Weyl ordering prescription for noncommuting operators, 
leading to the midpoint discretization in (|TT|). Our calculations were checked 



on the convergence with respect to the value of iV and usually N = 30 was 
found to be sufficiently large to detect no iV-dependence within statistical 
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T 

Fig. 1. Logarithmic derivative L(T) as a function of the time T for e 2 = 2.5m and 
[i = 10~ 3 m. The solid line is a fit of the data to the form (8). 



errors. We projected G on the eigenstates of total and angular momentum 
|P, m) but due to the nonlocality in (O) this did not imply that the degrees of 
freedom associated with the symmetries (P, L 2 ) could be integrated out. The 
logarithmic derivative L(T) can be symbolically written as a statistical average 
(here we denote the set of integration variables (s, s, {zi}, {zi}) shorthandedly 
as Z): 

_ fVZ {K'[Z\ - V'[Z\) e -K[z]+v[z\ 



jVZ e -K[z]+v[z] 



(12) 



where the prime stands for analytical differentiation with respect to T. We 
performed Metropolis Monte Carlo calculations of this object by averaging 
K'[Z] — V'[Z] over an ensemble generated by exp (—K[Z] + V[Z]). In Fig. 1 we 
present results of calculations of L(T) for the case e 2 = 2.5m and // = 10~ 3 m. 
It is clear that L(T) approaches a constant for large T and a fit of the data 
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Fig. 2. Mass Eq of the ground state for SQED2+1 as a function of the photon mass 
fi. The lines are just drawn to guide the eye. Any logarithmic dependence of Eq on 
\i would show up in this hgure as a straight nonhorizontal line. 



to the form (|§|) allows one to get an accurate result (within 1%) for E and 
in principle even a good indication (within 5%) of the first excited state. 
Note that we are able to do calculations at arbitrary large T since spacetime 
is not discretized in the FSR approach. This is in sharp contrast with lattice 
calculations where spacetime is necessarily finite. Besides that, the accuracy of 
our results is high as compared to lattice calculations |T(J . In Fig. 2 we display 
the mass of the ground state obtained with this procedure, as a function 
of fi for two different values of the coupling constant e 2 . Clearly the limit 
fi — > is a smooth one and we are able to restore gauge invariance in this 
way. Hence, no logarithmic confinement is found in (2+1) dimensions as long 
as the mass of the charged particles is kept finite. From these studies we 
furthermore find that for each /1 there is a finite critical coupling constant 
where E = 0. For \ijm = 0.1 this critical coupling constant was determined 
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to be e 2 rit /m = 6.0 ± 0.5. At this point the vacuum becomes degenerate in 
view of these zero mass two-particle bound states. 

To get some insight on the absence of the logarithmic confinement in the 



relativistic case we may study the BSE [I2|, r3"| , |l6| in the ladder approximation 
for this theory. In the CM system where P to t = (0, y/s) it takes on the following 
form: 

2 

$(p) = 7^ J d 3 g V(p, q) G (q) $(q) (13) 



with 



V(p,q) = f + ^T^r 2 aIld = T~2 1 1 2\2 i 2 ( 14 ) 

(p - q) 1 + /i 2 (m 2 - js + q 2 ) 2 + sq§ 



After Wick rotation JT4J the momenta in (|13|) and ( p^) can be taken Euclidean. 
The solutions to ( [13] ) were constructed by expanding $, V and G on the basis 
of 0(3) spherical harmonics Yi m (Q 3 ) and truncating the resulting infinite set 
of coupled one-dimensional integral equations at a certain Z max . Truncation 
at Z max = 2 was sufficient to get very accurate results. Further details of this 
procedure will be published elsewhere. In Fig. 2 the result of the calculations 
for e 2 = 2.5m are indicated by the diamonds. As compared to the FSR result, 
the BSE solutions show much less binding. This is a general feature of all our 
calculations and it is particularly striking for strong couplings. We are led to 
conclude that the crossed ladder diagrams give a very significant contribution 
to the binding energy of the system considered here. The limit \x — > is seen 
to be smooth as well and we can safely remove the IR cut-off this way. 

One may address the question what mechanism is responsible for the smooth- 
ness of the limit /i — ► in the relativistic calculations. It is generally known 
that both the BSE and the FSR have the correct nonrelativistic limit and 
besides that, it is rather remarkable that such a global feature as this loga- 
rithmic divergence does not persist when one considers a relativistic theory. 
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These considerations raise the question how the nonrelativistic limit is being 
reached in this case. For this purpose we carried out a variational analysis of 
the BSE for a g02+i~theory. It is convenient to introduce the nonrelativistic 
coupling constant A = g 2 /4m 2 which enters the corresponding Schrodinger 
equation. We verified that relativistically this case behaves smooth as well for 
H — > 0, while it becomes identical to sQED 2 +i in the nonrelativistic limit. Our 
trial function used: 



Viai(M) = exp 



12 2 
2 a T 



exp 



(15) 



with variational parameters a and /3, is particularly suited for evaluating 
the various matrix elements analytically. Nonrelativistically we expect a ~ 
O {^J Eb/tti^ ^> (3 ~ O (Es/m). Applying the Rayleigh-Ritz variational prin- 
ciple to the energy functional s(a,/3) indeed yields this property up to mul- 
tiplicative logarithmic corrections. The nonrelativistic limit is obtained by 
letting m/ /i — > oo and at the same time A//i — > 0, while keeping their product 
( = mX/fi 2 constant. In this region the wavefunction becomes independent of 
the relative time t and as a result the Schrodinger predictions are obtained. 
We find in particular that for 

-log 2 (m/A) <C' < 1 (16) 
m 

the binding energy Eg = (m 2 — ^sj jm goes as ~ ^Alog£. Hence in this 
region the logarithmic divergence of the Schrodinger analysis is recovered. We 
see however that Eb only exhibits this logarithmic dependence on the photon 
mass /i as long as ( is kept at a fixed value and m —>■ oo so that the condition 
( |T6| ) is satisfied. Taking on the other hand m large but fixed and letting fi — > 0, 
we effectively put ( —>■ oo, thereby invalidating (|TE|). It can be shown that in 
this limit E B approaches a large but finite constant ~ A log(m/ A) independent 
of /i. This shows that the nonrelativistic limit is not uniform. For decreasing 
fj, there is a crossover point /i (/io ~ A log(m/ A) as can be inferred from (|HJ)) 
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where the relative time dependence in the wavefunction starts to play a role. 
A dynamical screening mass /xq is effectively generated, which is related to 
the nonvanishing of the relative time parameter (3 in (|i~5l) . It is interesting to 
note that although /io is vanishingly small on the scale of m, it is essentially 
of relativistic origin and no trace of it is left in the Schrodinger equation. 

It is known that accurate results can be obtained variationally, even with 



rather simple trial functions ||15|| . Over a wide range of coupling constants we 
indeed find that the above variational calculations yield within a few percent 
the exact BSE results. In conclusion, we have shown that the FSR is a very 
suitable nonperturbative method to extract in a reliable way the bound state 
energy. Furthermore, the IR cut-off in sQED 2 +i can safely be removed. No 
logarithmic confinement as /i — > is found in (2+1) dimension for a finite 
mass of the charged particles. This is due to the relative time effects occurring 
in a relativistic description. 
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